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Abstract. Grigori Kuzmin was a very gifted dynamicist and one of the tower- 
ing figures in the distinguished history of the Tartu Observatory. He obtained a 
number of important results in relative isolation which were later rediscovered 
in the West. This work laid the foundation for further advances in the theory of 
stellar systems in dynamical equilibrium, thereby substantially increasing our 
understanding of galaxy dynamics. 

Key words: stellar dynamics - galaxies: elliptical and lenticular, cD - galax- 
ies: kinematics and dynamics - galaxies: structure 

1. INTRODUCTION 

At the scientific symposium celebrating the 200th anniversary of the Tartu Ob- 
servatory, it is appropriate to consider the work of one of its leading astronomers, 
Grigori G. Kuzmin. Kuzmin was born in 1917, was a student at Tartu University 
from 1936 to 1940, and became an assistant at the Observatory in the same year. 
He received his PhD degree in 1952 and continued to work at the Observatory until 
after his retirement in 1982. He passed away in 1988, having made fundamental 
advances in the theory of galaxy dynamics. 

In this contribution we concentrate on Kuzmin's work on galaxy models with 
separable potentials, from spherical to axisymmetric to triaxial (sections 2-4). The 
stellar orbits in these galaxy models all have three exact integrals of motion which 
are quadratic functions of the velocity components, allowing the calculation of 
many properties by analytic means. The separable models not only provide insight 
into the internal orbital structure and formation of galaxies, but they are also 
essential for testing more sophisticated numerical modeling approaches (section 5). 
Our conclusions follow in section 6. 
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2. DYNAMICAL MODELS FOR GALAXIES 

In order to understand the dynamics of a collisionless stellar system such as the 
Milky Way, the goal is to find the phase-space distribution function f(x, v) > 
which provides the distribution of stars over all positions x and all velocities v. 
The mass-density distribution p(x) follows by integrating f(x, v) over the veloc- 
ities: p = J fd 3 v. The velocities derive from the gravitational potential &(x) 
via Newton's equation of motion dv/dt = —VQ(x). If the stellar system is 
self-consistent, then the density is related to the potential by Poisson's equation 
4irGp(x) = V 2 $(f). 

Traditionally two main approaches are taken to construct dynamical models. 
Either one chooses a form for / in order to find p, which requires making a sensible 
Ansatz for / based on some pre-conceived notion on the galaxy formation process, 
or one assumes p and attempts to find / by solving an integral equation. In both 
cases it is advantageous to use Jeans' (1915) theorem which states that / depends 
on the six phase-space variables x and v only through the single- valued (isolating) 
integrals of motion admitted by the gravitational potential of the system. Since in 
non-spherical systems there are generally at most three such integrals, application 
of Jeans' theorem reduces the number of variables in the problem of constructing 
a dynamical model. 

The second order velocity moments (viVj) of / (with i,j — 1,2,3) are related 
directly to the density and the potential by the equations of stellar hydrodynamics, 
also called the Jeans equations. They can often be solved without knowledge of /, 
but this generally requires additional assumptions about the shape of the velocity 
distribution. Moreover, such Jeans models are not guaranteed to correspond to a 
distribution function / that is everywhere non-negative, i.e., physical. 

It is useful to consider spherical models. In this case the Hamilton- Jacobi 
equation separates in standard spherical coordinates (r,6,4>), and the motion of 
each star is governed by four isolating integrals of motion: the energy E, and the 
three components L x , L y and L z of the angular momentum vector. All orbits are 
planar rosettes (with the maximum angular momentum circular orbit and the zero 
angular momentum radial orbit as extremes). The mass distribution is defined by 
specifying the density profile p(r), and the associated gravitational potential $(r) 
can be found by carrying out two single weighted integrations of p(r). There is only 
a single Jeans equation. Many self-consistent models have been constructed by 
analytic means. These include isotropic models with / = f(E) found via an Abel 
inversion given the density p(r) (Eddington 1916), circular orbit models in which 
only orbits with zero radial action are populated, the Osipkov (1979)/Merritt 
(1985) models in which / = F(E + ah 2 ) and many more general distribution 
functions / = f(E, L) with L the total angular momentum. 

Over the past century many papers were written on the construction of spher- 
ical models. Popular mass models are Henon's (1959) isochrone and the 7-models 
(e.g., Dehnen 19933), which include the Jaffe (1982) and Hernquist (1990) models. 
Many of these were, in fact, studied much earlier by Kuzmin and collaborators, in 
particular Veltmann and later Tenjes. They derived density profiles, Jeans solu- 
tions and distribution functions. The results were not well known in the Western 
literature, but were summarized in English in Kuzmin (1993). However, Kuzmin's 
main contribution lies elsewhere. 



^^Also derived by M. Franx in unpublished notes dated 1988. 
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3. KUZMIN'S PIONEERING WORK 

When Kuzmin started working on galaxy dynamics in Tartu in the early nine- 
teen forties, under very difficult circumstances, the field concentrated on the Milky 
Way, for which many basic properties were still poorly known. This included 
the mass distribution, the associated gravitational potential, the nature of the 
three-dimensional stellar orbits other than those provided by the epicyclic approx- 
imation for nearly-circular orbits, and the phase-space distribution function. He 
contributed to all these topics. 

It was known from classical mechanics that in a flattened axisymmetric system 
such as the Milky Way, the energy E and the angular momentum component L z 
parallel to the symmetry axis [in cylindrical coordinates (R, </>, z)], are integrals of 
motion. If the phase space distribution function / would depend only on E and L z 
(by the Jeans Theorem), then the second moments of the stellar motions observed, 
e.g., near the Sun, should obey the relations (v\) — {v 2 z ) and {vrv z ) = 0. But it 
was already clear from the observational data that (vjj) ^ {v 2 z ) and {vrv z ) ^ 0. 
This strongly suggested that the distribution function / of the Milky Way depends 
on a third parameter, or in other words that the Galactic potential must support 
a third integral of motion /3 . The nature of this integral was not well understood, 
and at the time only an approximate expression been derived for nearly circular 
orbits with small vertical harmonic oscillations (Lindblad 1933). 

Non-spherical gravitational potentials with three exact integrals of motion E, 
I2 and J3 had been known for half a century (Stackel 1890, Eddington 1915, 
Clark 1937), but had received relatively little attention. The integrals I2 and 
^3 are quadratic functions of the velocity components. Chandrasekhar (1942) 
assumed that / = j[E + a/2 + W3) for such systems (the Ellipsoidal Hypothesis), 
and attempted to derive the associated density p. He demonstrated that the 
only self-consistent models satisfying this assumed form for / are spherical, and 
hence inconsistent with flattened systems such as the Milky Way. Interest in 
such spherical models was revived by Osipkov (1979) and Merritt (1985), but 
Chandrasekhar 's result did not encourage further work on the flattened separable 
models. There was also little interest in pursuing the opposite route: to try to solve 
for / for densities p that have a separable potential. Van Albada (1952) (wrongly) 
demonstrated that oblate axisymmetric separable potentials are not associated 
with mass distributions that are everywhere non-negative, suggesting that their 
usefulness for stellar dynamics was severely limited. 

In the mean time, Kuzmin had in fact made very significant progress in this 
area. He had worked on mass models for our Galaxy built by adding non- 
homogenoeus ellipsoids, work later repeated by e.g., Burbidge, Burbidge & Prcn- 
dergast (1962), had demonstrated that only isolating integrals of motion should 
be used in Jeans' Theorem, well before Lynden-Bell's (1962b) paper on the sub- 
ject, and, most significantly, had derived many interesting results for axisymmetric 
models with an exact third integral of motion. The results were collected in his 
seminal PhD Thesis (Kuzmin 19523), an d many of these were published in the 
decade that followed. 

Kuzmin considered axisymmetric mass models with a gravitational potential 

A — V 



2 Translated into English by Tenjes and Einasto, including additions from 1969. 
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where (A, (j), v) are prolate spheroidal coordinates, with (A, v) the two roots of 

R 2 z 2 

~^- + ^— = h (2) 

t + a t + 7 

and the constants a and 7 such that —7 < v < —a < A and F(t) an arbitrary 
smooth function (r = A, v). Surfaces of constant A are spheroids, and surfaces of 
constant v arc hyperboloids of two sheets, respectively. In these potentials, the 
Hamilton- Jacobi equation separates in (A, <fi, u), resulting in three exact integrals 
of motion E, L z and 7 3 , which are quadratic in the velocities and can be written 
explicitly. Kuzmin showed that such potentials have useful associated densities, 
given by a simple formula, and that the density p(R, z) > if and only if p(0, z) > 
0, a result now known as Kuzmin's theorem (de Zeeuw 1985b). 

Kuzmin assumed a simple form for the density profile along the short axis, 
p(0, z), where z 2 = t + 7, which in the notation used here would be equivalent to: 

^(r + 7) 1 / 2 F'(r)cxr-"/ 2 . (3) 
dr 

He showed that the case n = 3 provides a fair approximation to the Milky Way 
potential as it was known at the time (without dark halo); the model is nearly 
spheroidal and is in fact a flattened generalization of Hcnon's isochrone. The 
density (and potential) can be written explicitly in terms of R and z, avoiding 
the need to use spheroidal coordinates. The case n = 4 corresponds to an exactly 
spheroidal model p{R,z) = po/(l + m 2 ) 2 with m 2 = R 2 /a 2 + z 2 /c 2 and a = 
—a 2 , 7 = — c 2 . In the limit of extreme flattening both the n = 3 and the n = 4 
models reduce to an innnitcsimally thin disk with surface density S(i?) = So/(l + 
R 2 /a 2 f/ 2 . This is the Kuzmin (1953) disk which was rediscovered by Toomre 
(1963). It has many remarkable properties, including the fact that the potential 
above the disk is equal to that of a point mass at (R = 0, z = —a) below the disk 
and the potential below the disk is that of a point mass at (R = 0, z = +a) above 
the disk. Kuzmin furthermore noted that the model with n = no is the weighted 
sum of models with n > no; this built on his pioneering 1943 work on construction 
of models by superposition of inhomogeneous spheroids. 

The stellar orbits in oblate separable models are bounded by the orthogonal 
confocal coordinate surfaces, and are now called short-axis tubes. Their properties 
are similar to those of the orbits in the Milky Way found numerically by Ollongrcn 
(1962) using Schmidt's (1956) mass model. Kuzmin & Kutuzov (1962) calculated 
the f(E,L z ) for the model n = 3. They realized that p(R, z) can be written ex- 
plicitly as p(R, $) without any reference to spheroidal coordinates, which allows 
computing f(E, L z ) via series expansions a la Frickc (1952). This work occurred in 
parallel with Lynden-Bell's (1962a) construction of self-consistent flattened Plum- 
mer models using integral transforms. But Kuzmin & Kutuzov apparently did not 
investigate three-integral distribution functions for oblate (and the related prolate) 
density distributions. This was left to Dejonghe & de Zeeuw (1988) who made full 
use of the many elegant properties of the Kuzmin-Kutuzov mass model, derived a 
closed form for f(E, L z ) (see also Batsleer & Dejonghe 1993) as well as a number of 
three-integral distribution functions, and computed the projected surface density 
and the observed kinematic properties. 
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4. TRIAXIAL GALAXY MODELS 

Kuzmin also generalized much of his work on oblate separable models to those 
with triaxial shapes. This work is largely unpublished except for a very concise 
four-page summary (Kuzmin 1973) in the proceedings of the Alma Ata conference 
in 1972, which was translated into English in Kuzmin (1987). The paper con- 
siders triaxial models with potentials separable in confocal ellipsoidal coordinates 
(A, jU, v) defined as the three roots for r of 

^ +^ + — = 1, (4) 



t + a t + f3 r + 7 

with (x, y, z) the usual Cartesian coordinates, and with constants a, j3, 7 such that 
—7 < v < —f3 < n < —a < A. For potentials of the form 

$ s= m + m + m (5) 

3 (A - fi)(X — v) {ji- v){[i — A) [y — \)(y — p) {) 

all orbits have three exact integrals of motion E, I 2 and I 3 , quadratic in the 
velocities. Kuzmin noted that the orbits can be divided into four main families, 
that the density associated to $s is non-negative, i.e., p(x,y,z) > if and only 
if p(0, 0,z) > 0, and that it can be expressed elegantly via a 'Kuzmin formula'. 
The set of models includes an exactly ellipsoidal model, also known as the 'perfect 
ellipsoid': 

p{x,y,z) = p /(l + m 2 ) 2 , with m 2 = x 2 /a 2 + y 2 jb 2 + z 2 /c 2 , (6) 

where 7 = — c 2 , j3 ~ — b 2 and a = —a 2 . The four major orbit families are referred 
to as box orbits, inner long-axis tube orbits, outer long-axis tube orbits, and 
short-axis tube orbits, and are identical to the four major orbit families found 
in Schwarzschild's (1979) numerical model for a stationary triaxial galaxy with a 
constant density core (de Zeeuw 1985a). 

All these results for triaxial models were obtained via an independent route 
by the first author in the period 1982-1985, following a suggestion by Lyndcn- 
Bcll to 'look into how well separable potentials could be made to fit a general 
potential in the finite density core of a galaxy' (de Zeeuw & Lynden-Bell 1985). 
This lead to the discovery of the perfect ellipsoid (cq. 6) in a fairly roundabout 
way — as the ellipsoidal model for which the series expansion of the potential fits 
a separable one to all orders — and to the subsequent rediscovery of the Kuzmin 
(1956) paper on oblate separable models. The latter inspired the independent 
derivation of the Kuzmin formula for the triaxial case, and was followed by the 
full derivation of the gravitational potential for oblate, prolate and triaxial models 
given as the sum of two weighted single integrations of the density profile along 
the short axis (de Zeeuw 1985b; de Zeeuw, Peletier & Franx 1986). It was Jaan 
Einasto who cheerfully pointed out to the first author in 1985 during the Princeton 
IAU Symposium on Dark Matter that 'Kuzmin had already done it all for triaxial 
models many years earlier'. This overstated the case a little, as the derivation of 
the potential given the short-axis density profile was genuinely new. But the trail 
led to the Alma Ata conference contribution already mentioned, and the decision 
to have it translated into English and included as an Appendix of the proceedings 
of IAU Symposium 127 (Kuzmin 1987). The text was carefully checked by Kuzmin 
himself about one year before he passed away. 
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As Kuzmin knew, each separable triaxial model can be written as the weighted 
integral of constituent perfect ellipsoids, with the weight function found by means 
of a Stieltjes transform. It follows that the projected surface density is the same 
weighted integral of the projections of the constituent ellipsoids, which are ellip- 
tic disks. This property led to a new and powerful method to find the three- 
dimensional gravitational potential of non-axisymmetric disks (Evans & de Zeeuw 
1992). Franx (1988) had already deduced general properties of the projections of 
the separable models. All these properties are in fact shared by a much larger 
set of models: each ellipsoid p = po/(l + rn 2 ) p with p — n or p — n/2 generates 
a similar family (de Zeeuw & Pfenniger 1988; Evans & de Zeeuw 1992), but the 
motion is separable only for the models generated by the p = 2 ellipsoids. 

The second moment tensor is diagonal for the separable models, and the di- 
agonal elements are related to p and $5 by the Jeans equations, which are three 
partial differential equations for three unknowns (Lynden-Bcll 1960). Kuzmin and 
others had derived special solutions for the axisymmetric and elliptic disk limiting 
cases, but the general equations for triaxial geometry resisted solution for more 
than forty years, until van de Ven et al. (2003) found it by classical methods. As 
for all Jeans solutions, they are not guaranteed to correspond to a physical model 
with / > 0, but do allow efficient calculation of the second velocity moments. 

It has also been possible to construct analytic self-consistent models, but the 
existence of more than one major orbit family, in fact four, implies that f(E, I2, 13) 
is not uniquely defined by yo(A,/z, u), and certain assumptions have to be made. 
Populating only the tube orbits with zero radial action (thin tubes) together with 
the box orbits, Hunter & de Zeeuw (1992) constructed models that can be con- 
sidered as the triaxial generalization of the circular orbit spheres. Van de Ven et 
al. (2008) gave a comprehensive description of the so-called Abel models in which 
/ = J2i fi(E + a ih + bih), w hh and bi constants for each component i. These 
had been studied earlier by Dejonghe & Laurent (1991) and Mathieu & Dejonghc 
(1999), and are triaxial generalizations of the spherical Osipkov-Merritt models. 

The above brief summary shows that through Kuzmin's work and the sub- 
sequent follow-up, the theory of stationary triaxial dynamical models is now as 
comprehensive as that for spheres. This also means that by matching these mod- 
els with observations of galaxies — such as their observed surface brightness to 
infer the intrinisic density, and their observed stellar kinematics to constrain the 
gravitational potential and internal velocity distribution — we can gain insight 
in the internal dynamics of a much larger range of non-spherical stellar systems, 
including early- type galaxies. 

5. EARLY-TYPE GALAXIES 

Measurements of the stellar kinematics of early- type galaxies with integral- field 
spectrographs such as SAUR0N (Bacon et al. 2001) show that they can be divided 
into slow rotators and fast rotators, different from their morphological classifica- 
tion in ellipticals and lenticulars (Emsellem et al. 2007, 2011). Most fast rotators, 
including lenticular as well as many elliptical galaxies, seem to be close to oblate 
axisymmetric disk- like systems (Cappellari et al. 2007). On the other hand, vari- 
ous slow-rotator elliptical galaxies show clear deviations from axisymmetry in their 
photometry such as isophotal twists, as well as in their kinematics such as kine- 
matic misalignment and kinematically decoupled components (see e.g. Emsellem 
et al. 2004, and Fig. 1), indicating that they are triaxial. 
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Fig. 1. Triaxial dynamical model of the elliptical galaxy NGC4365. Left: Two- 
dimensional stellar kinematic observations of NGC 4365 obtained with the integral-field 
spectrograph SAURON (top row) accurately fitted with our numerical implementation of 
Schwarzschild's orbit superposition method in triaxial geometry (bottom row). From left 
to right: mean line-of-sight velocity and dispersion, and higher-order Gauss-Hermite ve- 
locity moments, which quantify deviations from a Gaussian velocity profile. Right: The 
black contours drawn on the maps are the isophotes derived from Hubble Space Tele- 
scope imaging of NGC 4365 shown on the right-hand side with a magenta box indicating 
the extent of the SAURON observations. Whereas the imaging shows a smooth ellipti- 
cal luminosity distribution, the kinematics reveal a complex non-axisymmetric orbital 
distribution with an (apparent) kinematically decoupled core. 



The separable triaxial galaxy models described above are able to capture much 
of the rich internal dynamics of these slow-rotator galaxies. However, the presence 
of central density cups and super-massive black holes imply that there are no 
global integrals I2 and I3. The tube orbits are relatively unaffected, but the box 
orbits are replaced by a mix of so-called boxlets (higher-order resonant orbits) and 
chaotic orbits (e.g., Gerhard & Binney 1985). While separable models are very 
useful to provide insights as well as tests (e.g., van de Ven et al. 2008), to model 
these elliptical galaxies in detail, a more general technique such as Schwarzschild's 
(1979) numerical orbit superposition method is needed. 

Schwarzschild's method allows for an arbitrary gravitational potential, with 
possible contributions from a central black hole and a dark matter halo. The 
equations of motions are integrated numerically for a representative library of or- 
bits. Next, the orbital weights are determined for which the combined density and 
higher order velocity moments best fit simultaneously the observed surface bright- 
ness and (two-dimensional) kinematics. The resulting distribution of (positive) 
orbital weights represents the distribution function (Vandervoort 1984), which is, 
by construction, guaranteed to be non-negative everywhere. Our recent imple- 
mentation in triaxial geometry (van de Ven et al. 2008, van den Bosch et al. 2008, 
van den Bosch & van de Ven 2009), allows the construction of accurate models for 
slow-rotator galaxies and to establish their internal phase-space structure. 

One of these galaxies is NGC 4365, a giant E3 galaxy with a core that appears 
to rotate around the short axis (Surma & Bender 1995; Davies et al. 2001), almost 
perpendicular (~82°) to the main body rotation around the long-axis (Fig. 1). 
The customary interpretation has been that this core is kinematically distinct (a 
so-called KDC) and a remnant of a last major accretion ~12Gyr ago. However, 
triaxial Schwarzschild models fitted to the SAURON stellar kinematics (van den 
Bosch et al. 2008; Fig. 1), revealed that the KDC is part of the single triaxial 
structure: short axis tube orbits dominate throughout, but half of them counter- 
rotate leading to vanishing net rotation around the short-axis, except in the center 
due to a small in-balance leading to the appearance of a KDC. The intrinsic shape 
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changes from oblate triaxial in the central regions to prolate triaxial in the outer 
parts where long-axis tube orbits provide the net rotation seen in the observed 
velocity field. 

Triaxial Schwarzschild models for another dozen of these slow-rotator ellipticals 
reveal a similar picture of galaxies dominated by short axis tube orbits with can- 
celing rotation apart from a small in-balance in the center leading to the KDCs. In 
addition, these slow-rotator galaxies contain on average 20% box orbits consistent 
with their mildly triaxial shape (van den Bosch et al., in prep.). Several fast- 
rotator galaxies also contain KDCs, but typically a magnitude smaller in size, and 
they seem to be truly decoupled axisymmetric disk-like components with a younger 
stellar population (see also McDermid et al. 2006). The slow- rotator galaxies with 
these apparent KDCs, in harmony with their smooth age and metallicity gradi- 
ents, provide a true challenge for (cosmological) merger simulations: Dry mergers 
over- predict the fraction of box orbits (e.g., Jesseit et al. 2005) but invoking gas 
seems to require special conditions to make the merger remnants as round as the 
observed slow-rotator ellipticals (e.g., Hoffman et al. 2011, Bois et al. 2011). 

6. CONCLUDING REMARKS 

Kuzmin also contributed to other areas of galactic dynamics, and trained a 
number of students, amongst which Einasto, Kutuzov and Veltman. The rela- 
tive isolation in which he worked resulted in a delayed recognition of many of his 
achievements in the West. This was due in part to limited possibilities for publish- 
ing, information exchange and travel in the Soviet era, and in part to provincialism, 
as few in the West were able to read Russian. Kuzmin's seminal 1956 paper was 
published in the Russian Astronomical Journal in the last year before it started to 
be translated in English. His papers often had a short English summary attached, 
and he sent English synopses of his work to key astronomers, amongst which 
Lynden-Bell, Oort and others, but these were not widely distributed. Perek's 
(1962) review did help advertize his results, but even so, some of his work was 
independently rediscovered. Today's on-line availibility of preprints and journals, 
including new and old Baltic Astronomy publications, and the strongly improved 
possibilities for travel, should go a long way towards avoiding that genuine discov- 
eries are overlooked and repeated. 

A number of the ideas in Kuzmin's (1956) paper on oblate separable models 
were generalized to prolate and triaxial models and in turn led to further develop- 
ments, as described in sections 3 and 4. These contributed significantly to bringing 
the modeling of separable triaxial systems to the level reached for spherical sys- 
tems: a theory of mass models and their projections, full understanding of the 
orbits, general solution of the Jeans equations, and construction of anisotropic 
distribution functions. This effort is valuable for testing more sophisticated nu- 
merical triaxial models and for gaining insight into the dynamical structure of 
early-type galaxies. As such, Kuzmin's work contributed substantially to our un- 
derstanding of galaxy dynamics. He belongs in the list of very famous astronomers 
associated with Tartu Observatory. 
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